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We report the results of a density-functional study of the diffusion of Pt dimers on the (111) surface 
of Pt. The calculated activation energy of 0.37 eV is in exact agreement with the recent experiment 
of Kyuno et al. [Surf. Sci. 397, 191 (1998)]. Our calculations establish that the dimers are mobile 
at temperatures of interest for adatom diffusion, and thus contribute to mass transport. They 
also indicate that the diffusion path for dimers consists of a sequence of one-atom and (concerted) 
two-atom jumps. 

PACS numbers: 68.35.Fx, 71.15.Nc, 71. 15. Mb, 68.35. Md 



I. INTRODUCTION 

The diffusion of adatoms — or very small clusters of 
them — on clean, infinite, defect-free surfaces remains, 
in spite of its conceptual simplicity, a largely unresolved 
problem. Experimentally, direct measurements of the 
trajectories of individual particles are possible (using 
field-ion microscopy) for only a few elements (notably 
Ir and Pt)El in a narrow range of temperatures. In most 
cases, the diffusivities are determined indirectly, e.g., by 
inferring them from island growth measurements (see for 
instance Ref. ||). The situation is just as difficult theo- 
retically: though the diffusion constants can in principle 
be determined explicitly using, e.g., molecular-dynamics 
simulations, this is in practice extremely difficult because 
of computer and model-potential limitations. 

The (111) surface of Pt is of particular interest. The 
diffusivity of adatoms has been measured using differ- 
ent approaches, and aji|-activation barrier of 0.25-0.26 eV 
has been determined an This is in_disagreement with our 
first-principles result of 0.33 eV,B but within the error 
bar of a recent calculation by FeibelmanjD which gives a 
barrier of 0.29 eV.El Though the agreement between the- 
ory and experiment can be judged satisfactory, questions 
remain concerning the role played by dimers in the ki- 
netics of growth, and in particular the shape of islands 
as a function of temperature. 

We have carried out a series of density-functional- 
theory calculations in order to determine the barrier and 
mechanism for diffusion of the Pt dimer on Pt(lll). 
The activation energy we obtain — 0.37 eV — is in 
striking (and somewhat surprising) exact agreement with 
experiments Our calculations therefore indicate that the 
dimers are certainly mobile at temperatures of interest 
for adatom diffusion (and thus contribute to mass trans- 
port), as in fact can also be inferred from the measure- 
ments of Kyuno et alB Our study reveals that the diffu- 



sion path for dimers consists of a sequence of one-atom 
and (concerted) two-atom jumps. 

II. COMPUTATIONAL DETAILS 

As already noted above, the calculations reported 
here were performed within the framework of density.-, 
functional theory.a Our previous studies of this surfaced 
indicate that the local-density approximation (LDA)Ed 
to the exchange-and-correlation energy provides a better 
description of platinuia than the generalized-gradient ap- 
proximation (GGA)jB the LDA was therefore used for 
all calculations reported here. The ion cores were ap- 
proximated by pseudopotentials with 5d electrons treated 
as valence states. The pseudopotentials were gener- 
ated using the semi-relativistic scheme of Troullier and 
MartinsEj and expressed in the KleinmaruBylander form 
using the s component as the local one.E3~t3 The elec- 
tronic wave-functions were represented using a plane- 
wave basis set with kinetic energy up to 40 Ry. In order to 
improve convergence, the electronic states were occupied 
according to a Fermi distribution with ksT c \ — 0.1 eV 
and the total energies obtained by extrapolating to zero 
electronic temperature. For similar reasons, the calcula- 
tions were initiated using wave-functions obtained from 
the self-consistent solution of the Kohn-Sham Hamilto- 
nian in a mixed basis set composed of pseudo-atomic or- 
bitals and plane waves cut off at 4 Ry.E3 The minimiza- 
tion of the energy with respect to the electronic [degrees 
of freedom was done using an iterative procedure.tU After 
achieving electronic convergence, the atoms were moved 
according to a damped Newton dynamics until forces be- 
came less than 0.01 eV/A. In view of the high energy 
cutoff needed in the plane-wave expansion, it is impor- 
tant to keep the system size to a minimum. To do so, 
we used a slab geometry consisting of four 3x3 layers 
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(plus vacuum in the supercell approach) with the dimer 
adsorbed on one side. Only the top layer (plus the dimer) 
was allowed to relax. The k-space integration was per- 
formed using a 2 x 2 grid, the exact number of points 
depending on the actual symmetry of the configuration 
under consideration. 



III. RESULTS 
A. Binding and dissociation energies 

The binding energy of the dimer is given by ^binding = 

-Edimcr + -^clcan — 2£' a dsorption, where -Edimcr IS the total 

energy of the system including the dimer, E c \ can is the 
total energy of the system with a "clean" surface, and 
^adsorption is the adsorption energy of a single adatom. 

A single adatom on the (111) surface can sit either 
in a fee site or in a (stacking-fault) hep-like site. The 
latter lies.II. 17-0.21 eV above the former in the case of 
platinumBu i.e., is very unfavourable. For the equilib- 
rium state of the dimer, one therefore expects the two 
atoms to sit in nearest-neighbour fee sites, such as the 
fi f2 configuration in Fig. ^. The corresponding config- 
uration where the two atoms are in hep-like sites, such 
as hih.2, indeed lies 0.30 eV above fif2, we have veri- 
fied. In both the fif2 and hili2 configurations, the dimer 
bond length is close to a/y/2 (with a the lattice parame- 
ter), the nearest-neighbour distance on the (111) surface. 
Other possible configurations of the dimer have longer or 
shorter bond lengths and are therefore unlikely. We have 
examined the fif3 configuration, where the dimer atoms 
are second nearest neighbours (bond length = \fia/\/2) 
and found it to lie 0.56 eV above fif2- Likewise, the pos- 
sible fep-hep configurations (such as fihi, fih.2, and fili3 
- cf. Fig. [y) all lie substantially higher in energy than 
fif2, as we will see below. 

With the dimer in the fjf2 configuration, we obtain a 
binding energy of 0.48 eV, with an error that we estimate 
to be of the order of 0.05 eV. This is much larger than the 
diffusion barrier for a single adatom, 0.33 cV.cl The bar- 
rier towards dissociation is given, roughly, by the sum of 
the dimes binding energy and the diffusion barrier of the 
adatomta — 0.81 eV in the present case. We therefore 
expect the dimers to be stable, i.e., unlikely to dissoci- 
ate, in the temperature range in which the adatoms are 
mobile. More precisely, the first step towards dissoci- 
ation is most likely related to the fif2 — » fib? (rather 
than fili2) transition (from geometrical considerations: 
the fili3 distance is larger than the fili2 one; also, there 
are two recombination paths for fih.2, but only one for 
fihs). We find the barrier for this process to be about 
0.75 eV, remarkably close to the above estimate. 



B. Diffusion by successive jumps 

We consider first the possibility that diffusion of the 
dimer proceeds by successive jumps of its constituent 
atoms. Referring to Fig. |l] again, we will assume that 
atom B, on the f2 site, jumps first, followed by atom A, 
on the fi site. B can jump to either of the three adjacent 
hep-like sites, as indicated by arrows in Fig. § The fthi 
(i = 1, 2, 3) configurations of the dimer are metastable 
states since hep sites are not equilibrium sites on this sur- 
face. Evidently, these states will be relevant to diffusion 
only if they lie sufficiently low in energy above the equi- 
librium state — more specifically by an amount which is 
of the order of the barrier for single-atom diffusion (0.33 
eV). 

We have calculated the energies of the fih,; configura- 
tions of Fig. [j], and obtain, as measured with respect to 
the fif 2 state, AE = 0.87, 0.34, and 0.69 eV for fihi, 
fih 2 , and fih 3 , respectively. Thus, of the three possible 
intermediate configurations, only fili2 is probable on a 
timescale comparable to that for adatom diffusion (but of 
course cannot be excluded on longer timescales). We note 
that the barrier for the fif2 — * fib-2 process — 0.35 eV 
according to our calculations — is very close to the fih.2 
configurational energy; thus, the barrier for the fili2 — > 
fif2 process is vanishingly small (an unsignificant 0.01 
eV) so that fih 2 is very short-lived. 

If diffusion proceeds via a sequence of single-atom dis- 
placements, then there are at this point two possibilities: 
(i) Atom A may jump to either hi or I14. (ii) Atom B may 
jump to either £4 or back to f 2 . It is easy to sec that, in 
either case, the next jump (of either A or B) would lead 
to an improbable high-energy state of type fihi or fih3, 
which both lie substantially above the equilibrium fif2 
state. We can therefore only conclude that diffusion by 
successive single-atom jumps only brings about a local 
motion of the dimer, which is essentially trapped in a 
potential well (possessing multiple minima) out of which 
it cannot escape on a timescale appropriate to adatom 
diffusion. 

We investigate next the possibility that diffusion pro- 
ceeds via the concerted motion of the two atoms forming 
the dimer. 



C. Diffusion by concerted jumps 

In order to go from one equilibrium site (such as fif2) 
to another in a concerted manner, the dimer must go 
through a metastable near-neighbour hep-hep configura- 
tion, as evident from Fig. [I]. As already noted above, the 
near-neighbour hep-hep configuration lies 0.30 eV above 
the equilibrium fcc-fcc state; this is comparable to the 
single-adatom diffusion barrier, and therefore is a possi- 
ble candidate for the dimer diffusion path. 

As depicted in Fig. [jj there exists three possibilities 
for such a concerted jump, labeled 'cji', 'cj2\ and 'cjs'. 
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The paths l cj2' and 'CJ3' are equivalent by symmetry, 
but different from 'cji': in the transition state of the 
cji path, the two atoms forming the dimer sit on either 
side of a surface atom, which is not the case for cj2 and 
CJ3. This makes diffusion extremely difficult in the cji 
direction; our calculations predict, indeed, a barrier of 
approximately 0.8 eV. In contrast, for the cj2 (or cquiv- 
alently CJ3) process, we obtain a value of 0.37 eV. This 
is only 0.07 eV above the hep-hep configuration energy. 
Thus, once in this state, the dimer can easily jump to 
an adjacent equilibrium fcc-fcc configuration (which can 
be either the initial one or a new one) via the cj4 or ej5 
process indicated in Fig. [I]. 

The barrier for the cj2 process, 0.37 eV, is very close to 
that for adatom diffusion; though we have not calculated 
the corresponding prefactors (they are found experimen- 
tally to differ by roughly an order of magnitude) jj it is 
therefore certainly the case that both processes will con- 
tribute significantly to diffusion at temperatures of inter- 
est. Kyuno et al.B indeed, find the diffusivity of dimers 
at 150 K to be comparable to that of adatoms at 100 K, 
with a difference in activation barriers of 0.11 eV in fa- 
vor of adatoms, larger thaiuthat observed here (0.04 eV). 
Following Bogicevic et al.jtB the temperature at which an 
Arrhenius process becomes active can be estimated from 
To = (i5 y 4/fcs)/ln(i/o/r) where Ea is the activation en- 
ergy, Vq is the attempt-to-diffuse frequency (prefactor), 
and r is the actual frequency at which diffusion is taking 
place. For dimers, Ea = 0.37 eV and vq sa 1.6 x 10j 2 s -1 
(using the prefactor determined experimentally) ;Lj with 
r«l s , corresponding to an experimental deposition 
rate of 0.001-0.1 monolayer per second, one finds that 
dimers become active at approximately 150 K. 

It is important to note that if only cj2-type processes 
are possible, then (because the barrier for cji-type pro- 
cesses is comparatively much larger) diffusion would be 
constrained to a one-dimensional corridor, consisting of 
a sequence of jumps such as Cj2-cj4-cj2-cj4... Correlated 
jumps nee d not, however, be exclusive of other processes. 
In section IIIB, we argued that successive jumps alone 
would not lead to mass transport, which does not mean 
that they do not contribute. In fact, once in a hep-hep 
configuration, the dimer can find its way to another cor- 
ridor by the combination of single-atom moves hili2 — » 
fih.2 — > fif2- The energies of the h]^ and fih 2 configu- 
rations are comparable (0.30 vs 0.34 eV) and the barrier 
for the fili2 — > fif2 process is a smallish 0.01 eV; this 
diffusion path is therefore highly probable. The corre- 
sponding barrier for the concerted hili2 — > fif2 process 
is, we have just seen, 0.07 eV. 

The error bar on the above values is of a few hundrcths 
of an eV and we therefore cannot determine precisely 
which route will be the preferred one. It can however 
be safely concluded that either will lead to significant 
mass transport: the limiting factor for diffusion is the 
fif2 — > hili2 barrier of 0.37 eV, quite comparable to that 
for single atoms. 



IV. CONCLUDING REMARKS 

Direct measurements of the diffusion of Pt atoms and 
dimers on Pt(lll) using low-temperature field- ion mi- 
croscopy have been reported very recently by Kyuno et 
alB The activation energy for adatoms is found to be 
0.260 ± 0.003 eV while it is 0.37 ± 0.02 eV for dimers. 
The barrier for adatoms is in excellent agreement with 
previous experimental cstimatesJj'Q which seems to rule 
out the possibility that the activation barrier for adatoms 
as deduced from growth experiments is "contaminated" 
by contributions from dimers. 

For dimers, the measured barrier is in remarkable 
agreement with our theoretical estimate. Such a close 
agreement is probably to some extent fortuitous as our 
calculations are precise to no more than a few hundreths 
of an eV because of our neglect of dynamical and quan- 
tum effects, finite-size limitations, and the approximate 
character of the LDA. Nevertheless, the present study 
does establish that dimers are mobile at temperatures 
where single-atom diffusion is active and can therefore 
contribute to mass transport, albeit perhaps not in a very 
significant manner compared to adatoms. Our calcula- 
tions, further, indicate that the pathway for dimer dif- 
fusion consists of a sequence of one-atom and concerted 
two-atom jumps. Such a diffusion mechanism has been 
reported recently in the case of Al dimers on Al(lll) 
based on density-functional theory calculations (see Ref. 
cited in ID. 
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FIG. 1. The dimer diffusion processes studied in the 
present work; the symbols are explained in the text. 
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